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Abstract 

We investigate the temperature dependence of the instanton content of gluon 
fields and their contribution to quark correlation using quenched lattice QCD 
and the cooling method. We found a suppression of the topological suscep- 
tibility at finite temperature, agreeing with the PCAC expectation at low 
temperature and enhanced suppression for temperatures at and above the 
deconfinement transition. For temperatures up to about 334 MeV, the topo- 
logical charge correlation agrees well with a single instanton profile, though 
the size parameter seems to change across the phase transition. The screen- 
ing wave functions for the pion and the rho become slightly more compact 
at higher temperatures. Lattice cooling shows no contribution from the in- 
stantons to the screening wave functions at temperatures above the phase 
transition. 

PACS numbers: 12.38.Gc, 11.10. Wx, 12.38.Mh 
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I. INTRODUCTION 



Lattice QCD calculations JTJ have suggested the existence of a deconfinement phase 
transition at a critical temperature T c , at which quark and gluons are deconfmed and form a 
new phase of matter, the so-called quark-gluon plasma. This possibility stimulated numerous 
theoretical and phenomenological works on the subject, many of which focus on the possible 
experimental signatures of the quark-gluon plasma 0j. Whereas these works are extremely 
useful, in view of the experiments at AGS, CERN, and RHIC, theoretical studies of the 
properties of finite temperature QCD matter should also be important. Lattice gauge theory 
provides a useful tool to calculate various hadronic matrix elements from first principle, and 
it therefore enables us to study static finite temperature properties of QCD matter. 

A particularly important and interesting question regarding QCD matter at finite tem- 
perature is what happens to the hadrons. Most of the phenomenological works on quark- 
gluon plasma signatures to date are based on a picture of weakly interacting quarks and 
gluons at temperatures above T c . However, such a simple picture is not completely sup- 
ported by lattice calculations. For example, lattice calculations of the so-called screening 
wave functions indicated strong correlations among the quarks in hadronic channels [[3]. 
However, it is still an open question whether the screening wave function results have any 
direct bearing on the structure of hadrons at finite temperature. 

In lack of a method to calculate the finite-temperature hadron wave functions directly, 
we turn our attention to a class of gauge configurations known to have important effects on 
the structure and masses of light hadrons at zero temperature ||], namely, large, isolated 
instantons. The primary goal of this paper is to study the instanton content of QCD matter 
and their contribution to quark correlations at finite temperature. 

It is well known that instantons break the {7a(1) chiral symmetry, resulting in the famous 
axial anomalies. At high temperature, Debye screening gives rise to an electric mass for the 
gluons and thereby suppresses instanton amplitudes; the Ua(X) symmetry should therefore 
be restored. Whether the C/a(1) restoration at Ty occurs at the same temperature as the 
SU(Nf)A restoration at T c bears important consequences to the hadronic physics at finite 
temperature: if T v = T c , the phase transition is of first order U and there are drastic 
changes to hadrons at T c ; if, on the other hand, Ty > T c , the phase transition is of second 
order. The two scenarios give rise to very different hadron physics at temperatures around 
T c ||. One of our goals is to locate Ty via lattice calculation of the topological susceptibility, 
which is a measure of the instanton content, and thereby differentiate the two scenarios. 

There have been calculations of the temperature dependence of the topological suscep- 
tibility in pure SU(2) 0,§, pure SU(3) 0, as well as unquenched SU(3) with two light 
quark flavors [jTH]. However, these authors vary the temperature by adjusting the coupling 
constant, making it tricky to compare the susceptibility at different temperatures because 
of the different lattice length scales. At larger (3 finite size effects may plague the calcula- 
tions using small (e.g., 8 3 x 4) lattices. There is also some confusion on the behavior of the 
susceptibility in the SU(2) theory. Whereas Teper et al. observed a sudden suppression of 
the topological susceptibility at the deconfinement temperature for both SU(2) and SU(3) 
fields Di Giacomo et al. reported that the cooling method gave ambiguous results 

for SU(2) theory, and the susceptibility determined by their field theoretical method stays 
almost constant across the phase transition J8|. 
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In this paper, we present calculations of the temperature dependence of the instanton 
content of quenched SU(3) fields using the cooling method. We vary the temperature by 
changing the number of time slices and report results for N t = 4, 6, 8, 10, 12, 14, 16, N s = 16, 
and (3 = 6 lattices. We did not observe the ambiguity reported for the SU(2) calculation, 
and the topological susceptibility is very stable with respect to the number of cooling steps. 
We found a suppression of the topological susceptibility starting at temperatures below 
T c . The topological charge correlation function agrees well with the continuum instanton 
profile, although there seems to be a change in the instanton size parameter across the phase 
transition. For three temperatures around T c , we also calculate the spatially-propagated 
wave functions, or the so-called screening wave functions, for the pion and the rho with 
and without cooling. We found that while the uncooled screening wave functions are only 
mildly sensitive to temperature, in agreement with previous calculations || , the cooled wave 
functions for T > 200 MeV are indistinguishable from the free ones. We thus conclude that 
at such high temperatures, instantons do not contribute to the screening wave functions. 

We first describe in Sec. II the set of observables we calculated. A brief outline of the 
cooling method is presented in Sec. Ill, and our main results are summarized in Sec. IV. We 
then close with some discussions in Sec. V. 



II. OBSERVABLES 

The instanton content of the gauge fields can be monitored by the the topological charge 
density, which can be defined on the lattice as 

Re Tr [U aP {x n )U l5 {x n )] , (1) 

OZ7T 

where U a p is the product of the link variables around a plaquette in the a — (3 plane, after 
some cooling steps were applied to remove lattice artifacts associated with discretization . 
The topological susceptibility is then given as the fluctuations of the topological charge: 



^<(EQM) 2 >, (2) 



Xt N t N':> 

where (. . .) indicates configuration averaging, and N t and N s are the number of sites in the 
temporal and spatial direction. The lattice spacing a chang function of the lattice 

inverse coupling (3, which makes it tricky to compare Xt calculated with different /3, especially 
when (3 is not large enough to be in the asymptotic scaling regime. To study the temperature 
dependence of the topological susceptibility, we chose to change the temperature by varying 
the number of time slices 

but keeping fixed (3 and hence a. This way the uncertainty in a enters only as an overall 
constant factor which does not affect the shape of the Xt vs. T curve. 

In order to look into the details of the topological charge distribution, we also calculated 
the topological charge density correlation function |4j] 
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c Q (x) = (Y,Q(y)Q(* + y))/(Y<Q 2 (y)) ■ ( 4 ) 

y y 

One can compare this correlation function to a convolution of an isolated instanton topo- 
logical charge density 

Qp( x ) = -n ( 2 P 2 ) ' ( 5 ) 

where p is the size parameter. Zero temperature calculations using the cooling method Q 
show that after about 50 cooling steps, the gauge fields are dominated by large, isolated 
instantons, whose profiles agree well with that given by the continuum analytic expression 
Eq. ([]) with p w 0.3 fm. We will monitor the temperature dependence of the topological 
charge density correlation as given by Eq. ([|). 

At zero temperature, due to Euclidean symmetry the space-like and time-like plaquettes 
(P x , P t ) have equal expectation values. At finite temperature, however, the Euclidean sym- 
metry is broken, and an asymmetry between the space-like and time-like plaquettes develops 
at the phase transition temperature proportional to the entropy s of the system ||11|| : 



Ts = e + P = 4(3 [1 - 0.16675/ + 0(g 4 )\ {(P t ) - (P x )) , (6) 

where 

(Pt) = — s— Y -TrU Qi (x n ) , (7) 

and 

s z n,i<j 

Since the pressure p is continuous across the phase transition, the discontinuity of e + p 
is the latent heat Ae. Understanding the physics origin of the latent heat is important to 
developing models for finite temperature hadronic matter. For example, it was postulated in 
Ref. ]12| that instanton liquid reorganizes into molecules near T c and this process contributes 



to a jump in the energy density at T c . Here we calculate the latent heat from both cooled and 
uncooled configurations so as to monitor the contribution of the instantons to this quantity. 

To study the quark distribution in various hadronic channels, we calculate the Bethe- 
Salpeter amplitudes |3"1,|T3"1| : 

x + y 

r _ if A(x')dx' 

^ S (y) = / dx(Vt\d{x)P e * -? 5 u(x + y)\n) (9) 

x+y 

r _ i J A(x")dx" 

V* s (y)= dx{Q\d(x)Pe * 7 2 u(f + y)\n) . (10) 

Here, \fl) is the vacuum state, and P indicates path-ordering of the gauge fields A put in 
to make the amplitudes gauge invariant. We have used the natural choice of a straight line 
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joining the two quarks at x and x + y . When T is finite, the above amplitudes are ex- 
tracted at large z instead of the usual Euclidean time filtering. This results in a projection 
into the lowest momentum state, and therefore one cannot interpret the resulting ampli- 
tudes as a wavefunction. However, some insights into hadronic structure can be obtained 
by comparing the QCD calculations to model results 0,0,0. In particular, one can com- 
pare the amplitudes for free quarks with those calculated at various temperatures, to check 
whether at some temperature, hadronic matter behaves like weakly interacting quark-gluon 
gas. Such a comparison has been done in two recent papers ||, with the surprising result 
that the Bethe-Salpeter amplitudes seems to be quite insensitive to changes in temperature 
up to T = 1.5T C . In this paper, we extend the previous calculations to higher temperature 
T = 2T C , and we also calculate the amplitudes from cooled configurations so that we can 
correlate any changes in the amplitudes with the instanton content of the gauge fields. 



III. COOLING 

The method of lattice cooling has been extensively studied and discussed in the literature 
||. The method can be applied to finite temperature configurations without modifications. 
We used the same procedure as in Ref. [f|]. 

Briefly speaking, the cooling method smooths out the gauge fields locally, by minimizing 
the action density. After some cooling steps the short-range fluctuations in the gauge fields 
are suppressed, leaving the long-range "bumps" more or less unchanged. In particular, large 
isolated instantons, which are protected by topology, dominate the gauge fields after many 
cooling steps, while the smaller instantons annihilate with nearby anti-instantons. The topo- 
logical charge approaches a plateau corresponding to an integral number of instantons when 
plotted against the number of cooling steps. The topological susceptibility then becomes 
constant, since the annihilation of instanton pairs does not change its value. Only then 
can one extract an unambiguous value for the susceptibility using the cooling method. In 
Fig. 1 we show the cooling history of (Q 2 ) for several N t . In contrast to Ref. ||, we have no 
problem extracting (Q 2 ) and therefore the susceptibility even for T > T c . 
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FIG. 1. The expectation value of topological charge squared vs. number of cooling steps for N t — 4, 
6, 8, 10, and 16. The quenched configurations are generated using f3 = 6 and N s = 16. In all cases, 
the values of (Q 2 ) become unchanged after only about 15 cooling steps. Same behavior was found 
for Nt = 12 and 14 as well, though not shown here for clarity. The error bars reflect statistical 
uncertainties in the configuration averaging only. 



After many cooling steps, the instantons are well separated in the gauge fields, and the 
dilute gas approximation should become valid. The probability distributions for topological 
charge is then given by a convolution of Poissonian distributions ||: 

i-(Q)=e- m «g^, (11) 

and the dilute gas value of (Q 2 ) is then given by 

(Q 2 > = 2m . (12) 

In the case of an ideal Poissonian process, m would be equal to the number of instantons 
(anti- instantons). A comparison of the distribution of topological charge at various Nt with 
the dilute gas approximation (dashed lines) is given in Fig. 2. Their agreement gives a 
nontrivial check that our gauge fields are dominated by dilute instanton ensembles after 
our cooling procedure, and the topological susceptibility we extract from Fig. 1 reflects the 
fluctuations in the topological charge in the ensembles. 
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FIG. 2. Histograms showing the distribution of topological charge in the configurations for N t — 4, 
6, 8, and 10. The dashed lines show the distribution given by the dilute gas approximation (see 
Sec. Ill), using the (Q 2 ) read from Fig. 1. Notice the clustering of data around integral values of 
the topological charge. 



IV. RESULTS 

For the generation of the gauge field configurations we used the Metropolis algorithm. We 
studied a lattice with spatial volume V s = 16 3 . Following Eq. (|3]), we varied the temperature 
by varying the number of time slices using values N t = 4, 6, 8, 10, 12, 14, and 16, respectively. 
The coupling constant is chosen to be (3 = 6.0 which corresponds to a physical lattice spacing 
a ~ 0.1 fm. Thus the calculation covers a temperature range between 125 and 500 MeV. For 
this value of (3 the critical temperature is T c ~ 250 MeV. We generated 100 configurations 
each for N t = 4, 6, 8, 10, 12, and 40 configurations for N t = 14 and 16, which were separated 
by 300 lattice sweeps. The configurations were then cooled for 40 steps using the Cabibbo- 
Marinari heat bath algorithm in the limit of infinite coupling strength b = oo as discussed 
in Ref. 01]. We verified the stability of our results by varying the number of cooling steps. 
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FIG. 3. Topological susceptibility as a function of temperature. Lattice data, shown with error 
bars, are compared with: 1) the PCAC expectation assuming soft pion gas as the heat bath, taking 
c = —1/6 and = 186 MeV (see Eq. 15), indicated by the solid line; 2) the Pisarski-Yaffe formula 
(Eq. 16, using p = 0.26 fm), derived assuming Debye-screening mechanism of instanton suppression, 
indicated by the dotdashed line; 3) the Pisarski-Yaffe formula rescaled to match with the value of 
the topological susceptibility at T c = 250 MeV, shown in dashed line. The phenomenological value 
is also plotted at zero temperature. 



The numerical values for the topological susceptibility are shown in Fig. 3. The plot also 
shows the approximate phenomenological value for Xt at zero temperature as estimated by 



Witten and Veneziano 



Xt 



f 2 

J IT 

4N~ f 



-2m 2 K + m 2 v ) ~ (180 MeV) 4 



(13) 



As can be seen, the phenomenological value and the numerical result for the V = 16 4 lattice 
agree quite well. The temperature behavior of Xt exhibits a gentle decrease for T < T c , 
turning into a relatively sharp decay of the susceptibility around the critical temperature 
T c . At low temperature, T <C T c , one can model the heat bath as a soft pion gas and apply 
PCAC methods to obtain 



Xt (T) = X t(T = 0) (l + cT 2 /F* 



(14) 



with —1/6 < c < 1/6. A very mild temperature dependence is thus predicted for low 
temperature. We show in Fig. 3 the maximum suppression of the topological susceptibility 
(c = —1/6) by the solid line, which agrees quite well with our data up to T w 200 MeV, 
though our statistics is not good enough to narrow down the range of a. The agreement 
is probably accidental though, since our quenched configurations should not reproduce the 
physics of a thermal pion gas. The instanton density is suppressed by a factor of one half 
when T = T c , with some stronger suppression mechanism operative at temperatures between 
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T = 200 MeV and T c . A rather sharp decrease of the topological susceptibility continues 
above the phase transition. At temperature T ~ 334 MeV topological effects are strongly 
suppressed; Xt is practically zero at T ~ 500 MeV. The dotdashed curve in Fig. 3 shows 



the temperature dependence for Xt as predicted by Pisarski and Yaffe 0] calculating the 
Debye-screening suppression of large scale instantons: 



XtiT) = Xt (T = 0) 1 + A73 2 exp 



-2A - 18a I + 7A 



(15) 



for quenched SU(3) fields, with A = npT, a = 0.01289764, 7 = 0.15858, and p = 0.26 fm 
is the size parameter of instantons we obtained by fitting our N t = 6 data (see discussion 
below), which differs slightly from a previous lattice calculation at zero temperature J3J. The 
perturbative result in Eq. [15] is supposed to be valid only at very high temperature. In the 
temperature range we are considering, the perturbative formula gives too large a suppression. 
The dashed curve in Fig. 3 was computed by rescaling the Pisarski- Yaffe expression to match 
with the susceptibility at T c = 250 MeV, which then seems to be in agreement with the data, 
though more data is needed between T c and 2T C for a more meaningful comparison. The 
phenomenological consequences of the restoration of [/^(l), such as changes in the 77' mass, 
should be an exciting subject for experimental works at AGS, CERN, and RHIC. 
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FIG. 4. Topological charge correlation function, calculated for Nt — 6, 8, 10, 12, 14, and 16. The 
solid lines show the correlation function calculated with the continuum instanton charge distribution 
assuming a size parameter p = 3.3a. The dashed line shown at Nt = 6 is the continuum instanton 
correlation function with p = 2.6a. 



In Fig. 4 we show the correlation function Cq(x), Eq. (|J), with x taken in the spatial 
direction. Six different temperatures are shown. For comparison the plots include the 
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correlation function calculated with an isolated instanton with size parameter p = 3.3a ~ 
0.33 fm, which fits the data reasonably well for all temperatures up to T c , at N t = 8. At 
T ps 334 MeV, or N t = 6, the instanton profile with p = 2.6a ~ 0.26 fm gives a better fit 
to the data. In Fig. 5, we show the cooling dependence of the topological charge correlation 
function for N t = 6. The shape of the correlation is very stable between 20 and 100 cooling 
steps, and the smaller size parameter p = 0.26 fm is clearly favored by the data over the 
larger one (p = 0.33 fm) preferred at lower temperatures. As can be seen, the temperature 
dependence is quite weak below T c , but our data suggests a rather sudden change in p across 
the phase transition. This may signal major reorganization in the gluon fields above T c . For 
N t = 4, the topological charge correlation function changes continuously as a function of 
cooling steps, and therefore we could not extract a meaningful shape at that temperature. 
The correlation function becomes flat at small N t when x is taken in the temporal direction 
because of the periodic boundary condition. 
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FIG. 5. Topological charge correlation function for Nt — 6, at different cooling steps. The solid and 
dashed lines are the same as in Fig. 4. 
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FIG. 6. Plaquette expectation values as a function of temperature for (a) cooled and (b) uncooled 
configurations. The space- like plaquette expectation values (P x ) (squares) are compared to the 
time-like ones (Pt) (diamonds) and the overall averages (P) (crosses). The lines are drawn to guide 
the eyes only. 



We monitor the contribution of the large, isolated instantons to the latent heat by cal- 
culating the splittings of the spatial and temporal plaquettes for both cooled and uncooled 
configurations. The numerical values of the plaquettes are plotted in Fig. 6 (a) and (b), 
where a clear splitting between the spatial and temporal plaquettes, AP = (P t ) — (P x ), can 
be seen for temperatures above the phase transition. The splitting AP, which is propor- 
tional to the latent heat, is clearly a function of cooling, and it approaches zero in the limit 
of infinite cooling steps. An estimate of the fraction of latent heat contributed by instantons, 
/, is given by 

/AP coo l cc l / -i r>\ 

» XB • ( l6 ) 

i-*-* uncooled 

In Fig. 7 we show / as a function of cooling steps for both N t — 4 and 6. In both cases, / 
drops very quickly in the first 40 cooling steps, showing that most of the splitting is due to 
short range fluctuations that got eliminated by cooling. After about 50 cooling steps, the 
N t = 4 splitting levels off to a plateau, where its value represents an upper bound to the 
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contribution of the instantons at T = 500 MeV. The splitting at N t = 6 has not reached a 
plateau yet even at 100 cooling steps; we quote just the fraction / at 100th cooling step as 
the upper bound at T = 334 MeV. Our results for / are 0.15 ± 0.04 at T = 334 MeV and 
0.056 ± 0.004 at T = 500 MeV. 
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FIG. 7. Plaquette splittings, AP C = (P t ) 
of cooling steps. Nt — 4 data are shown as crosses, and Nt — 6 data as diamonds. Values at 100th 
cooling step represent upper bounds to the fraction of latent heat contributed by the instantons. 



Finally, we computed the screening wave functions for the ir and p mesons as explained in 
Sec. m We used a conjugate gradient procedure to calculate the Greens functions for Wilson 
fermions with k = 0.1541, which corresponds to a quark mass m ~ 100 MeV. The average 
of a set of ten wave functions for N t = 10,6,4 corresponding to T = 200,334,500 MeV 
are shown in Figs. 8 and 9. As already observed in |3j the screening wave functions change 
only slightly as a function of temperature across T c . The spatial extension of the wave func- 
tions shrinks as temperature increases, in agreement with theoretical arguments discussed 
in Refs. Here the main idea is that in the deconfmed phase, the quark mass for 

propagation along the spatial axis is essentially given by the lowest Matsubara frequency 
m ~ ttT. Strong correlations resulting from unscreened magnetic interactions generate a 
bound state between quarks propagating along a spatial axis even in the deconfmed phase. 
The effective quark mass increases at high temperatures, thus decreasing the characteristic 
size of the state. However, the weak dependence of the wave functions is surprising in view 
of the factor of 2.5 increase in the effective quark mass from N t = 10 to 4. It is probable that 
the present temperature range is not high enough for the above picture based on dimensional 
reduction to be accurate. 
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FIG. 8. Screening wave functions of the ir meson, for N t — 4 (squares), 6 (diamonds), and 10 
(crosses), corresponding to T = 500, 334, 200 MeV, respectively, wave functions calculated with 
cooled configurations and symmetrized to take into account of the periodic boundary conditions are 
also shown for N t = 6 (circles) and 10 (small squares), which agree well with free wave functions 
(dashed and solid lines). 




FIG. 9. Screening wave functions of the p meson, see Fig. 8. 
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As in the case of zero temperature, cooling effectively removed short range fluctuations in 
the gauge fields responsible for quark confinement as well as the Coulomb force. As a result, 
the wave functions calculated from the cooled configurations are distinctly broader, with 
the slope at short distances strongly reduced (§]. To study the contribution of instantons 
to the screening wave functions at such high temperatures, we compare the cooled wave 
functions with the uncooled ones as well as the free wave functions. As shown in Figs. 8 
and 9, for both ir and p at T = 200 and 500 MeV, after symmetrized to take into account 
the periodic boundary conditions, the cooled wave functions are indistinguishable from the 
free ones. Our results imply that the large scale instantons do not play an important role 
in the screening wave functions in this range of the temperature. This is possibly a result 
of both the suppression of instanton density and the increase of the effective quark mass 
by the Matsubara frequency. Instead, some non-trivial interactions are operative, which 
correlate the effective massive quarks. It is important to understand this mechanism in 
order to develop an accurate phenomenology of hadronic physics at temperatures around 
T c , just where the current and future heavy-ion experiments probe. 

V. SUMMARY 

We discussed the content and the role of instantons in quenched QCD matter at finite 
temperature. Specifically, we calculated the density of instantons and found that it is sup- 
pressed as temperature increases. While the mild dependence of Xt can be understood by a 
PCAC method assuming a soft pion gas as the heat bath, the high temperature (T > T c ) 
behavior deviates significantly from the perturbative formula expressing the Debye-screening 
suppression of instantons. Our results indicate that the restoration of Ua{^) symmetry is a 
continuous process starting below T c ; at T c , about half of the axial anomalies are already 
removed. Some non-perturbative mechanism enhances the rate of L r y i(l)-restoration in the 
temperature range T c < T < 2T C . The topological charge correlation function agrees well 
with the continuum instanton profile, though our data suggests a rapid change in the size 
parameter from p = 0.33 fm at and below 250 MeV to 0.26 fm at 334 MeV. By comparing 
the entropy in the cooled and uncooled gauge fields, we obtain upper bounds — of about 5% 
at 2T C and 15% at 1.5T C — on the instanton contribution to the latent heat. The screening 
wave functions of 7r and p mesons show a weak temperature dependence for temperatures as 
high as 2T C . This points to strong correlations between quarks even in the deconfined phase. 
In contrast to the zero temperature situation, the screening wave functions at temperatures 
above T « 200 MeV are not dominated by the instantons. 

It will be interesting to compare these results with unquenched calculations, as it is quite 
possible that the reorganization of the gauge fields near the phase transition is sensitive to 
the presence of dynamical quarks. Much more theoretical work is needed to illuminate the 
complex and drastic changes of hadronic physics near the phase transition. 
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the U.S. Department of Energy, Grant No. DE-FG02-87ER40365, at Indiana University. 
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